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Abstract. We analyze the phase diagram of two-flavor quark matter under neutron star constraints for a 
nonlocal covariant quark model within the mean field approximation. Applications to cold compact stars 
are discussed. 

PACS. 12.38.Mh, 24.85. +p, 26.60. +c, 97.60.-s 



1 Introduction 

The characteristics of the QCD phase diagram is presently 
an important open subject of research in particle physics. 
In particular, the behavior of strongly interacting mat- 
ter in the region of low temperatures and large baryon 
densities could be tested against observational constraints 
from neutron stars PQ. Since in this thermodynamical re- 
gion one finds strong difficulties when trying to perform 
lattice QCD calculations, the pictures emerging from dif- 
ferent effective models of strong interactions deserve great 
interest from the theoretical point of view. Here we focus 
on a two-flavor chiral quark model which includes covari- 
ant nonlocal four-fermion interactions, motivated by an 
effective one-gluon exchange (OGE) picture. Nonlocality 
arises naturally in the context of several successful ap- 
proaches to low energy quark dynamics, such as the in- 
stanton liquid model [2] and the Schwinger-Dyson resum- 
mation techniques 3], and it is also a well-known feature 
of lattice QCD gj. 



2 Formalism 

The Euclidean action for the nonlocal model considered 
here, in the case of two light flavors and anti-triplet di- 
quark interactions, is given by 



Se = / d^x < tp(x) (— ift + m) ip(x) 



Here m is the current quark mass, which is assumed to 
be equal for u and d quarks, and Jm,d are mesonic and 
diquark nonlocal currents. The nonlocality is introduced 
here in a covariant way, through a separable interaction 
arising from an effective OGE picture. In this way, we have 

3u( x ) = J <^ z 9(z) i>(x + |) r f tp(x - |) , 

3d{x) = J d 4 z g(z) 4> c {x + |) i lb T 2 X a ip(x - |) , (2) 

where 4> c (x) — 7274 V" 7 \x) , // = (1,^751"), and r and A a , 
with a = 2,5,7, stand for Pauli and Gell-Mann matrices 
acting on flavor and color spaces, respectively. The func- 
tion g(z) is a form factor that characterizes the nonlocal 
interaction. 

The effective action in Eq. {l} might arise via Fierz 
rearrangement from some underlying more fundamental 
interactions, and is understood to be used — at the mean 
field level — in the Hartree approximation. In general, the 
ratio of coupling constants H/ G would be determined by 
these microscopic couplings; for example, OGE interac- 
tions lead to H/G = 0.75. Since the precise derivation of 
the effective couplings from QCD is not known, here we 
will leave H/G as a free parameter. 

Standard bosonization of the theory leads, in the mean 
field approximation, to a thermodynamical potential per 
unit volume given by 



(1) 



qmfa = 



\A\ 



2G 2H 



OO 

E 

n— — 00 



d 3 P 

(2tt) 3 



In det 



s- 



T 



, (3) 



2 



D. Gomez Dumm et al: Phase diagram of neutron star quark matter in nonlocal chiral models 



where the inverse propagator S~ 1 (a, A) is a 48 x 48 ma- 
trix in Dirac, flavor, color and Nambu-Gorkov spaces (its 
explicit form is given in Ref. [5]). Here a and A stand 
for the mean field values of scalar meson and diquark 
fields. Owing to the nonlocality, they come together with 
momentum-dependent form factors. The values of a and 
A can be obtained from the coupled gap equations 



df2 A 



dA 



, 



da 



= 



(4) 



In general one has to consider a different chemical 
potential /i/ c for each quark flavor / and color c. How- 
ever, when the system is in chemical equilibrium, not all 
fif c are independent. In our case, it can be seen [5] that 
they can be written in terms of only three quantities: the 
baryon chemical potential /i b , the quark electric chemical 
potential hq and the color chemical potential hs . Defining 
H = Mb/3, the corresponding relations read 



Hgb = H + Qq fXQ -2fl 8 /3 



(5) 



where, q = u,d, and Q q are quark electric charges. 

In the core of neutron stars, in addition to quark mat- 
ter we have electrons. The latter can be thermodynami- 
cally treated as a free Fermi gas, and their contribution 
has to be added to the grand canonical thermodynamical 
potential. Moreover, quark matter and electrons have to 
be in [3- equilibrium. Thus, assuming that antineutrinos 
escape from the stellar core, we must have 



Hdc Hue 



-HQ = He 



(6) 



If we now require the system to be electric and color 
charge neutral, the number of independent chemical po- 
tentials reduces further: He and hs are fixed by the condi- 
tions of vanishing electric and color densities. In this way, 
for each value of T and h we should find the values of A, 
<7, He and H8 that solve Eqs. (j4|), supplemented by the /?- 
equilibrium and electric and charge neutrality conditions. 



3 Numerical results 

According to previous analyses carried out within nonlo- 
cal scenarios [6], the results are not expected to show a 
strong qualitative dependence on the shape of the form 
factor. Thus we will consider here a simple, well-behaved 
Gaussian function 



g(p 2 ) = exp( ~p 2 /A 2 ) 



(7) 



where A is a free parameter, playing the role of an ultra- 
violet cut-off. 

For definitcness, for the input parameters we choose 
m = 5.12 MeV, A = 827 MeV and GA 2 = 18.78. These 
values are fixed so as to reproduce the empirical values 
for the pion mass and decay constant f„, and lead to 
a phenomenologically reasonable value for the chiral con- 
densate [7] at vanishing T and Hb, namely (OjgglO) 1 / 3 = 



—250 MeV. The only remaining free parameter is the cou- 
pling strength H in the scalar diquark channel. We choose 
here values for H/G in the range from 0.5 to 1, i.e. around 
the Fierz value H/G = 0.75 discussed above. 
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Fig. 1. Phase diagrams for the nonlocal OGE-based model, 
for different values of the ratio H/G, under neutron star con- 
straints. 



Our results for the phase diagrams are shown in Fig. [TJ 
where we plot the phase transition curves on T — /i dia- 
grams for different ratios H/G. In the graphs we show 
the regions corresponding to different phases, as well as 
the position of triple points (3P) and end points (EP). 
Besides the region of low T and Hi m which the chiral 
symmetry is broken (CSB), one finds normal quark mat- 
ter (NQM) and two-flavor superconducting (2SC) phases. 
Between CSB and NQM phases one has first order and 
crossover transitions — represented by solid and dotted 
lines, respectively — , whereas between NQM and 2SC re- 
gions, in all cases, we find a second order phase transition 
— dashed lines in the diagrams of Fig. [TJ Close to this 
NQM/2SC phase border, the dashed-dotted lines in the 
graphs delimit a band that corresponds to the so-called 
gapless 2SC (g2SC) phase. In this region, in addition to 
the two gapless modes corresponding to the unpaired blue 
quarks, the presence of flavor asymmetric chemical poten- 
tials gives rise to another two gapless fermionic quasipar- 
ticles. For the range of parameters considered here, how- 
ever, the g2SC region is too narrow to lead to sizeable 
effects. Finally, for intermediate values of H/G we find 
a region in which a certain volume fraction of the quark 
matter undergoes a transition to the 2SC phase coexist- 
ing with the remaining NQM phase. This is a 2SC-NQM 
mixed phase in which the system realizes the constraint 
of electric neutrality globally: the coexisting phases have 
opposite electric charges which neutralize each other, at a 
common equilibrium pressure. 

It can be seen that the 2SC phase region becomes 
larger when the ratio H/G is increased. This is not surpris- 
ing, since H is the effective coupling governing the quark- 
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quark interaction that gives rise to the pairing. As a gen- 
eral conclusion, it can be stated that even under neutron 
star constraints, provided the ratio H/G is not too low, 
the nonlocal scheme favors the existence of color supercon- 
ducting phases at low temperatures and moderate chem- 
ical potentials (we do not find 2SC only for H/G — 0.5). 
This is in contrast to the situation in the NJL model [8], 
where the existence of a 2SC phase turns out to be rather 
dependent on the input parameters. Our results are also 
qualitatively different from those obtained in the case of 
noncovariant nonlocal models [5], where above the chiral 
phase transition the NQM phase is preferable for values 
of the coupling ratio H/G < 0.75, and a color supercon- 
ducting phase can be found only for H/G w 1. It is now 
desirable to extend the present studies to other pairing 
channels relevant for neutron star cooling, such as spin-1 
pairing, where at present only the NJL [lOj and nonco- 
variant nonlocal models [llj have been considered. 



4 Application to cold compact stars 

One of the most important present applications of the mi- 
croscopic approaches to quark matter is to study whether 
such a state of matter can exist in the interior of cold 
compact stars. Let us briefly discuss our ongoing inves- 
tigations on this issue. Unfortunately, a consistent rela- 
tivistic approach to the quark-hadron phase transition, 
where hadrons appear as bound states of quarks, has not 
been developed up to now. Thus, we apply a so-called two- 
phase description, in which the nuclear matter phase is 
described within the relativistic Dirac-Brueckner-Hartree- 
Fock (DBHF) approach (see e.g. Ref. [12]) and the tran- 
sition to a quark matter phase is obtained by a Maxwell 
construction. From the curves presented in the previous 
section, it can be seen that in our OGE-inspired nonlo- 
cal model one finds a relatively low value of the critical 
density at T=0, hence some extra repulsion is needed in 
order to obtain a more realistic value. We have found that 
this can be achieved by including some small additional 
interaction in the omega vector meson channel, which does 
not affect in general the qualitative features of the phase 
diagrams discussed in the previous section. Given the cor- 
responding equation of state, the mass and structure of 
spherical, nonrotating stars is obtained by solving the 
Tolman-Oppenheimer-Volkov equations. Our preliminary 
results confirm the findings of Ref. [13] . i.e. that compact 
stars with quark matter cores are consistent with modern 
observations. Moreover, for the nonlocal quark models de- 
scribed here, this is found to be possible for values of H/ G 
closer to the standard, OGE motivated ratio H/G = 0.75. 



5 Conclusions 

We have studied the phase diagram of two-flavor quark 
matter under neutron star constraints for a nonlocal, co- 
variant quark model within the mean field approximation. 
The form of the nonlocal coupling has been motivated by a 



separable approximation of the OGE interaction. We have 
considered a nonlocal form factor of a Gaussian shape, 
and the model parameters (current quark mass m, cou- 
pling strength G, UV cutoff A) have been fixed so as to 
obtain adequate values for the pion mass, the pion decay 
constant and the chiral condensate at vanishing T and fis ■ 
After the numerical evaluation of the gap equations 
at finite temperature and chemical potential, considering 
different values for the coupling strength in the scalar di- 
quark channel, we have found that different low-tempera- 
ture quark matter phases can occur at intermediate densi- 
ties: normal quark matter (NQM), pure superconducting 
(2SC) quark matter and mixed 2SC-NQM phases. A band 
of gapless 2SC phase appears at the border of the super- 
conducting region, but this occurs in general at nonzero 
temperatures and should not represent a robust feature 
for compact star applications. Finally, in the context of 
the nonlocal theory discussed here, we have obtained pre- 
liminary results showing that compact stars with quark 
matter cores turn out to be consistent with modern ob- 
servations. 
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